00 

o 

O 
(N 

u 
< 



Scalability of quantum computation with addressable optical lattices 

T. R. Beals 
Department of Physics, University of California, Berkeley, California 94720, USA 

J. Valc(3 and K. B. Whalejjt] 

Department of Chemistry and Pitzer Center for Theoretical Chemistry, 

University of California, Berkeley, California 94720, USA 

(Dated: April 11, 2008) 

We make a detailed analysis of error mechanisms, gate fidelity, and scalability of proposals for 
quantum computation with neutral atoms in addressable (large lattice constant) optical lattices. 
We have identified possible limits to the size of quantum computations, arising in 3D optical lattices 
from current limitations on the ability to perform single qubit gates in parallel and in 2D lattices 
from constraints on laser power. Our results suggest that 3D arrays as large as 100 x 100 x 100 sites 
(i.e., ~ lO'' qubits) may be achievable, provided two-qubit gates can be performed with sufficiently 
high precision and degree of parallelizability. Parallelizability of long range interaction-based two- 
qubit gates is qualitatively compared to that of collisional gates. Different methods of performing 
single qubit gates are compared, and a lower bound of 1 x 10~^ is determined on the error rate for 
the error mechanisms affecting ^^^Cs in a blue-detuned lattice with Raman transition-based single 
qubit gates, given reasonable limits on experimental parameters. 
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I. INTRODUCTION 

Neutral atoms trapped in optical lattices constitute 
a promising system for quantum information process- 
ing [H [2]. Single qubit operations and qubit readout 
have already been demonstrated [3], albeit in a non- 
scalable system, and a number of two-qubit gates have 
been proposed [H [5]. Addressable optical lattices — in 
which the lattice spacing is large enough that individual 
lattice sites can be targeted by a laser [B] — offer an en- 
vironment that can be scaled to thousands of qubits in 
a three-dimensional (3D) array |7]. Preparation, load- 
ing, and imaging of an addressable optical lattice have 
recently been demonstrated [H]. 

Achieving large-scale fault-tolerant quantum computa- 
tion requires single and two-qubit gates with extremely 
high fidelity, as well as the ability to perform many gates 
in parallel. Current estimates of the fault tolerance er- 
ror threshold range from 10~^ to 10~^ for conventional 
quantum error correction, depending on the difficulty of 
communication between physically distant qubits and the 
ability to prepare certain states "offline" in a reliable 
manner [5J [TUl HH HU US]- More radical error correc- 
tion schemes [Ml [151 US! Hi] may offer better thresholds, 
but at the cost of high overhead. For the optical lattice 
system that is the focus of this work, the most relevant 
estimates [12] suggest a fault tolerance threshold on the 
order of 10~^. 

While much research has focused on schemes for re- 
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alization of qubits and performing quantum operations 
in a wide variety of experimental systems, the de- 
tailed physics and scalability of specific architectures 
have to date received less attention. Some work has 
been done on system-level analysis of architectural is- 
sues [m [m EHl El] . However, to provide real numbers 
to questions such as how large a system of qubits may be 
made and how many quantum operations can be made 
on this, it is necessary to undertake detailed analysis of 
the behavior of the proposed qubits in situ. The different 
physics involved in different implementations pose a va- 
riety of challenges, some of which may be system-specific 
while others, such as achieving gate fidelities with fault 
tolerant threshold values are quite generic. For exam- 
ple, the need to use thermal ensembles rather than pure 
states — as is usual in gas state and solid state proposals — 
has presented a major challenge for liquid state NMR 
because of limitations imposed by small thermal polar- 
ization and difficulties of initialization [251 1131 IMl US]. 

For qubits defined in internal states of trapped neu- 
tral atoms or trapped ions, both the qubit interac- 
tions and their environmental decoherence mechanisms 
are well understood [2B]. This enables architectural is- 
sues to be examined with full microscopic analysis of all 
physical features of such qubits. For ion traps, analy- 
sis of the limiting features of trapped ion physics led 
to the proposal of a modular, multiplexed trap archi- 
tecture to allow scale up from a few ( < 10) ions to 
many, possibly thousands of individually addressable 
ions [27l [28i 20, and progress demonstrating components 
of such an architecture for small numbers of ions is now 
underway [SD] [311 [Ml [331 [31] • For neutral atoms trapped 
in optical lattices, up to 250 individual atomic qubits 
have been trapped and imaged in an optical lattice sys- 
tem that can be readily scaled up to include thousands 



of atoms and whose spacing is large enough to allow in- 
dividual addressability [8] . Large arrays of sublattice ad- 
dressable trapped atoms have also been made [351 ES]- 
Trapped neutral atoms and ions are somewhat com- 
plementary; individual addressing and quantum gates 
are more straightforward to implement with ions and 
have already been demonstrated experimentally for small 
numbers of ions (e.g., forming entangled states of < 8 
ions [371 135] and implementing algorithms with two [33] 
and three qubits [40]) but achieving scaleup to even hun- 
dreds of ions still presents a serious technical challenge. 
In contrast, although scaleup to lattices containing hun- 
dreds of neutral atoms has been demonstrated in an ad- 
dressable system [8] and gates have been demonstrated 
in non-addressable systems [35] and in addressable dipole 
traps [HJ 135], they have not yet been accomplished in 
lattices containing individually addressable atoms. Very 
recently, techniques for addressing individual sites in a 
lattice with small spacing have been proposed [43 , 44J . 

In contrast to the situation for atomic and ionic qubits, 
progress in scaleup of solid state realizations of qubits — 
such as Josephson junctions or "^^P in Si — from pairs to 
many qubits is in a far more rudimentary state. The kind 
of detailed analysis that is required to develop specific 
physical devices cannot be made yet, since the underly- 
ing microscopic physics of the qubits when in situ is not 
currently well enough understood, although significant 
progress is being made |45| . Furthermore, while solid 
state systems are often generically referred to as 'scal- 
able' because of the ability to fabricate large scale solid 
state devices, the individual elements or qubits are not as 
reproducible as gas phase qubits due to the complexities 
and variability of their surroundings 135]. Nevertheless, 
architectural studies are beginning to be made for these 
systems [47ll48] . 

In this work, we go a step further in analysis of scaleup 
for trapped neutral atoms, undertaking a detailed phys- 
ical investigation of the effects of single qubit errors and 
other imperfections that limit the scalability of neutral 
atom quantum computation in an individually address- 
able optical lattice. We consider different candidate sin- 
gle qubit gates and their sensitivity to various sources 
of experimental error. We then compare this to calcula- 
tions of the threshold rate for fault tolerant computation 
under the appropriate conditions, in order to estimate 
how large a quantum computation may be made within 
current technological constraints and possible near-term 
improvements. To our knowledge, this is the first such de- 
tailed physical estimation of a practical limit on physical 
scaleup for any proposed experimental implementation 
of pure state quantum computation. We hope that this 
detailed analysis for trapped neutral atoms will spur sim- 
ilar analyses for other physical implementations once the 
relevant microscopic physics is better understood. Given 
that no experimental system will have unlimited scalabil- 
ity, such physical analysis of technical limits to scalable 
systems of functioning qubits within current technology 
is an essential complement to theoretical algorithmic scal- 



ing characteristics derived from complexity theory. 

We restrict our analysis here to addressable optical lat- 
tices. While we do not explicitly consider the alternative 
lattices with global addressing that are also being stud- 
ied experimentally [49], we shall make some comments at 
the end of this paper on relative benefits that these other 
lattices might offer. The analysis in this paper employs a 
combination of perturbation theory and numerical tech- 
niques such as the pseudo-spectral method with a Cheby- 
chev decomposition of the Schrodinger propagator [50J to 
quantify the effects of both memory and gate errors de- 
riving from all known sources for trapped atoms. In some 
respects our calculations complement and extend those of 
Saffman and Walker ^1^ for *^Rb atoms in dipole traps. 
However, that work did not address the issue of scalabil- 
ity that is analyzed here after the various error rates have 
been quantified. When a specific choice of parameters is 
necessary, we consider here -'^^■^Cs atoms in an address- 
able optical lattice [8] of lattice constant a = 5 ^m, with 
depth of Ul = 200 fiK. The lattice is ortho-rhombic in 
geometry, and is created by blue-detuned beams at 800 
nm, with an intrapair angle of '^ 9° for each of the three 
pairs of beams. For the one-dimensional case, the lattice 
potential is given by: 
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Field-insensitive sub-levels of the 6s ^6*1/2 hyperfine 
ground-state manifold are chosen as the qubit basis: 
|0) = |F = 3,mi. = 0) and |1) ee |F = 4,mi. = 0). 
The auxiliary levels |2) = |F = 4, mi? = l) and |3) = 
|i^ = 3, rap = 1) are also involved in the single qubit gate 
presented here (see Fig. I2I. The procedure for loading 
and initializing the lattice is described in detail in Weiss 
et al [7], and Vala et al [SI]. We assume that atoms 
are cooled to the vibrational ground state, e.g., with 3D 
Raman sideband cooling [53J [54] . 



This rest of this paper is organized as follows. Sec- 
tion jll A| contains an analysis of the error mechanisms 
due to the lattice itself, such as scattering processes and 



loss of atoms from the lattice. Section II B analyzes the 



single qubit gate proposal based on a Raman two-photon 
process and Section [II C| analyzes microwave pulse-based 
single qubit gates. We consider the effects of off-resonant 
transitions of non-target atoms, scattering, heating of 
target atoms, and addressing beam targeting and inten- 
sity errors for both types of single qubit gate. Section 
jIIDj contains a brief comparative discussion of the scal- 
ability of different classes of two-qubit gates. Section 
jlllj provides an analysis of the results from the previous 
sections and their collective implications for the scala- 
bility of quantum computation in an addressable optical 
lattice. In section IV we summarize our conclusions, dis- 



cuss some possible ways to bypass the limitations iden- 
tified here and identify some useful applications within 
the constraints established here. 



II. ERROR MECHANISMS 

Fault tolerance thresholds are sometimes expressed in 
terms of a "unified" error rate comprising both storage 
and gate errors, but are often also written in terms of 
separate gate and storage error rates (and occasionally 
also preparation and read-out error rates). In order for 
error correction to be able to keep up with storage errors, 
a practical quantum computer must be able to perform 
many gates in parallel. Consequently, estimates of error 
threshold values have typically assumed that gates can 
be performed on arbitrarily many qubits in parallel |13j . 
Storage errors can then be considered on a similar footing 
to gate errors that occur with a frequency given by mul- 
tiplying the storage error rate by the typical gate time, 
Ti to obtain an effective 'error per gate time' that can 
be combined with the error per gate (EPG) in analysis 
of overall error rates. 

Fault-tolerance threshold theorems assume maximal 
parallelizability |13j , implying that all or nearly all qubits 
can be addressed simultaneously {ua — N, where N is 
the total number of physical qubits and n^ is the num- 
ber that may be simultaneously addressed) . In most pro- 
posed schemes for quantum computing this is extremely 
hard, except for the trivial case where one desires to per- 
form the exact same gate on all atoms simultaneously. 
Parallelizability thus constitutes an important figure of 
merit, since if ha grows slower than N, the effective stor- 
age error rate will eventually exceed the capacity of any 
error correction protocol. If only a fraction ua/N of 
qubits can be addressed, the effective EPG for a storage 
error will be approximately equal to the storage error rate 
mu ltiplied by the ratio NTi/ha- As we show in Section 
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UA is on the order of N^/"^ in 3D lattices, while in 2D 
lattices ua can be on the order of N. Ti is on the order of 
tens of microseconds for the microwave pulse-based sin- 
gle qubit gate, whereas with sufficient laser power, it can 
be nanoseconds or less for the Raman single qubit gate. 
In the remainder of this section we derive expressions 
for the EPG for various decoherence mechanisms encoun- 
tered by atoms trapped in an addressable optical lattice, 
using ^'^■^Cs as a specific example where necessary. 
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FIG. 1: (Color online) The solid blue line shows the Raman 
scattering cross-section o-(aj_t) for a ^^^Cs atom in a blue- 
detuned optical lattice, as a function of lattice light wave- 
length, Xl- The dashed green line shows the frequency- 
dependent polarizability a{uiL) at the frequency ujl = 
2ttc/Xl- Both a{u>L) and a{uiL) are calculated for e+i polar- 
ized light interacting with the initial state |_F = 3,mF ~ 0), 
with results for the \F — 4, mp = 0) state or opposite polar- 
ization light being similar. 



which they do not, are sources of storage errors. Fortu- 
nately, the decohering effects of Rayleigh scattering can 
be partially suppressed with pulse sequences [55]. For 
Raman scattering, no such method exists, and so we fo- 
cus our analysis here on this form of scattering. The 
effective storage EPG due to Raman scattering is given 
by ■:^'^[-, where Ti is the gate time and T/h the scat- 
tering rate. Calculating the relative transition strengths 
(see Appendix the details), we find that roughly half of 
the errors induced by Raman scattering will be bit-flip 
errors, while the rest will be leakage to non-qubit states. 
The Raman cross section u{liJl) (where oj^ is the lattice 
light frequency) can be calculated as described in the 
Appendix, and is shown in Fig. [l] as a function of the 
lattice light wavelength. The scattering rate is related to 
the cross section by 



A. Optical lattice-induced storage errors 

We assume that the lattice has already been prepared, 
and that each lattice site is initially occupied by exactly 
one atom in the motional |0) state. A detailed description 
of a procedure to achieve this perfectly loaded lattice is 
contained in Vala et al [52] . 



1. Photon Scattering 

Both Raman scattering, in which the initial and fi- 
nal states of the atom differ, and Rayleigh scattering, in 
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where E"^ is average over an atomic spatial distribution 
of the peak electric field squared. We note that the op- 
tical lattice potential depth is given by Ul{ujl,Eq) — 

-^ \a{cjL)\ J91 . The polarizability is shown as a func- 
tion of wavelength in Fig. IT] For an atom in the ground 
state in a red-detuned lattice E^ ~ Eq, whereas in a blue- 
detuned lattice, W = 2i^¥^El, where cj^ = f yJW^Jm 
is the characteristic trapping frequency. We can then 
calculate the Raman scattering rate for the blue-detuned 



and red-detuned cases: 
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Using Eq. (Sal, we see that, for •^■'Cs in a blue 



detuned optical lattice with the reference parameters 
given in Section III we obtain a Raman scattering rate of 
T/h = 2.2 X 10 s~^, and thus an effective EPG value 
of (2.2x10-^ s-i)^^. 



2. Qubit loss and leakage 

Qubit loss errors are particularly serious, in that they 
can not be automatically corrected by error correcting 
codes. When an atom is lost from the optical lattice, 
or leaks into a non-qubit state, it is necessary to first 
detect the error before it can be corrected. The lost atom 
must be replaced before standard erasure error correcting 
codes [56] can be applied to correct the error. Detecting 
qubit loss requires that we have a method of detecting 
the presence of an atom at a given lattice site without 
disturbing its state. 

Preskill [5T| identified a simple circuit for performing 
such loss detection measurements. The circuit requires 
an ancilla in a known state, two applications of a CNOT 
or CPHASE gate, a similar number of single-qubit gates, 
and a measurement of the ancilla. This measurement 
could fail by giving an incorrect result (false positives or 
false negatives), or by disturbing the state of the target 
atom. The latter type of error could be corrected by 
standard error correcting codes, while the former could 
be minimized by repeating the measurement as necessary. 
Another possibility for detecting qubit loss involves the 
use of a cavity QED system [58] . 

The need for having certain atoms in the lattice re- 
served for use as ancillas for this scheme could be 
avoided by transporting an extra-lattice ancilla atom 
where needed through the use of optical tweezers |59j . 
If an atom loss was detected, this ancilla would already 
be on hand to serve as a replacement. A drawback of this 



approach is that performing such operations in parallel 
would require many sets of optical tweezers. In the case 
of most leakage errors, parallelizable methods exist for 
detecting "leaked" atoms and returning them to a qubit 
state. 

Fortunately, qubit loss rates are very low, with storage 
times as long as 25 s already reported [F . Collisions with 
background gas atoms are the primary cause of loss, and 
so it appears that storage times can be increased further 
through improved vacuum systems. It is also possible 
that a method may be found for performing loss detec- 
tion measurements in parallel, which, when coupled with 
a means for replacing lost atoms, would allow qubit loss 
errors to be handled by standard error correction tech- 
niques. Consequently, qubit loss is not likely to be the 
dominant source of errors in the near future, and we will 
not consider it further in this paper. 



B. Raman- based single qubit gates 

Two-photon Raman transitions present an attractive 
option for single qubit gates because of the associ- 
ated speed of qubit manipulation. Raman-based sin- 
gle qubit rotations have recently been experimentally 
demonstrated on a time scale less than 100 ns for a sin- 
gle ^^Rb atom trapped in an optical dipole trap [¥T| . 
A theoretical analysis of factors contributing to gate im- 
perfections for a single ^^Rb atom concluded that gate 
fidelities of ~ 10~^ are possible [ST]. We analyze here 
the error mechanisms arising during Raman gates imple- 
mented for ^■^■^Cs atomic qubits in a blue-detuned optical 
lattice. 

We consider a Raman process in which the QSi/2{F — 
3) <-> 6P1/2 transition is driven with strength fii by e+i- 
polarized light at a detuning of Ai, and the QP1/2 *-* 
QSi/2{F = 4) transition is driven with strength il2 by 
e+i-polarized light at a detuning A2. 

In the general case of a Raman-based single qubit 



gate with two-photon detuning A = Ai — A2 



m. 



|f}i|^)/2(Ai + A2), Rabi frequency flR ~ r2ir2^/(2Ai), 
and pulse duration t, we have an effective off-resonance 
Rabi frequency ft' ~ -y/jilflp -I- A^, and the rotation is 
approximately described by the following matrix f51j: 
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Unless otherwise noted, we assume zero two-photon de- 
tuning, i.e., A = 0. For the specific resonance case 



Ai ~ A2, the rotation matrix is as follows: 

f cos(0/2) ze-^'f>sm{9/2) 
^ '^' ^ I ie'"^ sin(6l/2) cos(6l/2) 



with = \^R\t and = arg(riij). 

It is necessary to define a metric for fidelity of rotation 
operations. We consider a qubit in an arbitrary initial 
state V' undergoing a rotation R(6,(f>), and compare it 
to a TT pulse, Roin, 0). The fidelity is then given by the 
following relation ( 5b I : 
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where the overline represents an average over initial 
states and the corresponding error is given by P = I — F. 
(Note that this definition differs from 51J which consid- 
ered the fidelity of a tt/2 pulse on a specific initial state.) 



1. Neighbor atom errors 

In the case of a single qubit gate performed with two 
orthogonal Raman lasers, an atom that is adjacent to 
the target atom and that is on the axis of one of the 
two lasers will experience a small undesired rotation. 
The effective Rabi frequency for this non-target atom is 
ri'^ — QjiI{a,a)/Io, where I{r,z) is given by Eq. (19 1. 
From Eq. (5a I, we can determine the fidelity error in 



the desired identity operations for the four neighboring 
non-target atoms as: 
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where Xr ~ 894 nm is the Raman laser wavelength 



t{SUi — SUo)/h. For a tt pulse, t 

e+i polarized light we have S(f) = ,, , ^^ 

the range 50 GHz < Ai/27r < 5000 GHz, the ratio 

SUi-SUo 



= tt/I^rI, so for 
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We now wish to calculate the variance var((50) = 
{S(f>'^) — (Scj))'^ due to atomic motion and spatial variation 
in the Raman beam intensity. Since, for typical param- 
eters, the Rayleigh length of the Raman beam will be 
much larger than the beam waist {zq 3> wq), we need 
consider only motion in the transverse direction. In the 
transverse direction at the beam waist, the intensity has 
the form /(r, 0) = Ioe~^^ ^^'> , where /q is the intensity at 
the center, wq the beam waist, and r the transverse dis- 
tance from the center. The atomic motional states can 
be approximated by eigenstates of the two-dimensional 
harmonic trapping potential obtained by parabolic ex- 
pansion of the transverse potential at the minima of the 
lattice potential. 

For an atom in the resulting two-dimensional har- 
monic oscillator eigenstate \nx, Uy), we calculate the vari- 
ance of the phase shift using the fourth-order Taylor ex- 
pansion of the Gaussian beam intensity [55", /(r, 0) ~ 
/o(l — 2r^/w§ -I- 2r'*/wQ). This results in the variance 
var„^ ,„ ((50) 
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(In studying the temperature dependence of this effect, 
the reader may find it helpful to make the approximation 
kBT/2 ~ huj{n^ + 1/2) ~ huj{ny + 1/2).) From the ex- 
pression for fidelity of a tt pulse, Eq. ( 5b I , we see that the 
expected error probability will be P,a = |var„^.^„^(50). 



2. Spontaneous emission 

In the limit where Ai ~ A2, the probability of spon- 
taneous emission during a tt pulse is Pj-s — 2\a i ' "^ti^re 
T = 34.9 ns is the natural lifetime of the 6P1/2 state. For 
|r2/{| << |Ai|, the detuning Ai is related to the Rabi 
frequency and Raman laser intensity / by 
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therefore the probability of spontaneous emission is 

7: \nR\ 



Pr: 



2/t(8.3- 1012 W-1 m2 s-2)- 
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3. Raman beam AC Stark Shifts 

The difference in AC Stark Shift between the logi- 
cal |0) and |1) states gives rise to a phase shift S4> = 



4. Atomic motion 

In addition to the effects discussed above, atomic mo- 
tion will introduce noise through variation in the effective 
pulse area, |f2i?|i, and variation in the two-photon detun- 
ing, A. The former effect is simply a result of atomic 
motion across the Gaussian profile of the Raman beams, 
and has a similar form to the result calculated in the pre- 
vious section, Eq. (l9|. For a it/2 pulse, we obtain the 
following result for the variance: 
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This variation in pulse area will then result in an error 

Prpa = ivar„^,„^(|rifl;|i) for a tt/2 gate. 

Doppler shifts of the Raman beams will cause varia- 
tion in the two-photon detuning A. Unlike the isolated 
two-site dipole trap situation considered by Saffman and 
Walker [ST], our system involves a 3D lattice and thus 



does not allow for a convenient first-order Doppler-free 
Raman laser configuration. We thus expect to see signif- 
icant variation in the two-photon detuning due to atomic 
motion-induced Doppler shifts, as described by the fol- 
lowing relation: 
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From the general expression for the rotation matrix, 
Eq. (HI), we determine that this variation will result 



in a fidelity error Pr„ 
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5. Polarization effects 

The Raman beams used to perform the single-qubit 
gate have a Gaussian profile. This means that, even at 
the beam waist, the beam will have a small component of 
polarization other than the desired e+i, according to |51j : 
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where z/j — 7rwQ/2 is the Rayleigh length. This extrane- 
ous polarization can result in leakage errors by causing 
transitions to states outside the computational basis. To 
estimate the probability of such errors, we determine the 
relative magnitude of the second term above, which corre- 
sponds to undesired polarization "seen" by a target atom 
in a vibrational state jn^;, Uy). Eq. ( 13 1 suggests that this 



can be estimated by the ratio of the spatial extent of the 
atom to the Rayleigh length, i.e., by 
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where lUt — - \/'^UL/rn is the characteristic trapping fre- 
quency in the harmonic approximation. There are four 
possible non-basis states into which the qubit could leak: 
\F = 3,mF = ±1) , \F = 4,m,p — ±1). Since the matrix 
elements for the unwanted transitions are comparable to 
those for the desired transitions |60j , and the correspond- 
ing Rabi frequencies are smaller by a factor of IPolrei, 
the probability of leakage into any particular state for 
a TT gate can be approximated by sin (|7'o|rci7i'/2) — 
(|'7^o|rci7r/2) . The total leakage probability in the mo- 
tional ground state, Ux = riy = 0, is then four times that 
quantity, i.e.: 
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6. Laser intensity noise and line-width 

Noise in the Raman lasers affects the fidelity of the 
gate. If the relative intensity fluctuation is SI /I, then 
by Taylor-expanding eqn. (5b) with 9 ^ 0o + S9, we see 



this will result in an initial state-averaged fidelity error of 
P = (l/6)(7r(5///)2 for a tt gate and P = (1/6) {^51 /if 
a 7r/2 gate. 

Even if the Raman lasers are actively stabilized, shot 
noise provides a lower bound on relative intensity fluctu- 
ations. The fluctuation due to shot noise [27] is given 

by 
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where t] is the quantum efficiency of the detector used 
in the stabilization circuit, ajj^ is the frequency of the 
Raman laser, and Pr is the laser power. If we assume 
r] = 0.5 and that our stabilization circuit reaches the 
lower bound, then the minimum fidelity error is 
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For detunings much smaller than the absolute optical fre- 
quency, we can use Wi^ ~ 27r • 3.5 • lO"^'' s~^ (i.e., the 
Raman transition frequency). For a Gaussian beam of 
power Pr/2, the intensity / at the waist is related to the 
power as / = Pr/ttwq, while t.^ = ir/flu, with fl^ given 
by (|7|. This results in the following estimate for laser 
intensity fluctuation-induced error: 
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C. Microwave-based single qubit gates 

Site-specific single-qubit gate operation in 3D lattices 
can also be achieved through the use of a far off-resonance 
addressing laser focused on a single lattice site, combined 
with pulsed global microwave fields [6j [52l [6TJ [62]. In 
order to address a single atom, a Gaussian beam with 
waist substantially smaller than the lattice spacing is 
used, which results in the target atom seeing a much 
greater field than any neighboring atom. The intensity 
of a Gaussian beam is described by 

I{r,z) = /o"^exp(-2rVu''(^)) (19) 

where wq is the beam waist, Iq the intensity at the center 
of the waist, zq = ttwq/X is the Rayleigh length, and 

w(z) = woa/1 + ^ the beam width as a function of the 

V ^0 

axial coordinate z. 

The addressing beam causes an AC Stark Shift A^^ of 
the various levels of the target atom. Here we consider 



the scheme for ^^^Cs that is outlined in Figure |2] By 
choosing the "magic wavelength", Am, for the addressing 
beam (for ^^SCs, Aa/ :^ 880 nm), the |2) and |3) auxil- 
iary levels receive AC Stark Shifts of equal magnitude but 
opposite sign, while the qubit levels |0) and |1) are unaf- 
fected. This allows the |0) ^ |2), |2) ^ |3), and |1) ^ |3) 
transitions to be driven by global microwave pulses that 
are resonant only for the target atom. Alternatively, a 
coUimated beam could be used to address an entire row 
of atoms (provided the row was not much longer than the 
Rayleigh length ttwq/Xm), allowing ua = N^^^ identical 
operations to be performed in parallel. However, since 
the beam waist wq must be smaller than the lattice spac- 
ing a, this implies that ua < tta/Am and consequently 
only a relatively small number of atoms can be addressed 
simultaneously using this method. We discuss this limi- 
tation further in the next subsection. 

We have developed a software package, quantum sim- 
ulation software (QSIMS) [63 , for simulating the quan- 
tum dynamics of one and two-qubit gates in this and 
other systems. Using QSIMS, we discretize the spa- 
tial wave function of the atom on a grid, with a sepa- 
rate grid representing each possible internal state of the 
atom. Quantum dynamics are simulated by applying the 
Schrodinger propagator, expanded in Chebychev poly- 
nomials [50 . The kinetic portion of the Hamiltonian 
is applied by means of a Fourier transform of the dis- 
cretized wavefunction from the position basis to the mo- 
mentum basis. Transitions between levels are treated 
with a dressed state approach ^64j. Although QSIMS is 
capable of simulating three spatial dimensions, the sym- 
metry of the system and the near-separability of the lat- 
tice potential make it reasonable in most cases to perform 
simulations in only one or two spatial dimensions. This 
results in a significant speed-up, since the run time of the 
simulations is 0{N\ogN) for a grid of N points. 

We simulate the microwave pulse-based single qubit tt 
gate with QSIMS, using the parameters a = 5 /im, Ul — 
200 ^jK, Aac = 0.2 MHz, wo ^ 1.2 fim, and Ti = 76 ^^s 
(with microwave pulse intensity fii = 41341 s""'^ cho- 
sen appropriately to achieve this gate time). The atom 
is assumed to be initially in the |0) qubit state and in 
the motional ground state, so the final state of the gate 
correspond to |1) and the motional ground state. We 
investigated various diff'erent versions of the gate with 
QSIMS and found that the best performance is achieved 
by simultaneously driving all three transitions |0) -s-^- |2), 
|1) ^^ |3), and |2) ^^ |3), with pulse intensity chosen 
such that the strength U2 of the first two transitions is 
•\/3/2rii, where fli is the strength of the |2) <-> |3) tran- 
sition. 



1. Off-resonant transitions 

Since the microwave pulses proposed to perform these 
single qubit operations are applied globally, i.e., to the 
entire lattice, there is a small probability that such a 
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FIG. 2: (Color online) Schematics of a single qubit flip 
operation. A focused far off-resonant addressing laser pro- 
duces an AC Stark Shift Aac of the |_F — 4,mF = 1) and 
\F = 3,mF = 1) levels in a single target atom. The qubit 
levels are the mp = levels, labeled |0) and |1). Levels |2) 
and 1 3) are auxiliary levels involved in the single qubit gate. 
(Most of the hyperfine sub-levels not involved in the gate are 
not shown here.) Three global microwave pulses tuned to the 
transitions |0) ^^ |2), |1) <-> |3), and |2) ^^ |3) drive these 
transitions in the target atom with transition strengths f22, 
fl2, and fii. Other atoms are not affected to first order (see 
below) since their transitions are not resonant. Our simu- 
lations indicate that best performance is achieved by simul- 
taneously driving all transitions, with the relative transition 
strengths chosen such that ^2 = \/3/2f2i. 



pulse will cause a given non-target atom to undergo a 
non-resonant transition. We can minimize this probabil- 
ity by carefully tuning the gate parameters such that the 
pulse ends with non-target atoms in a local minima of 
their Rabi cycles, and by making the detuning. A, large 
compared with the Rabi coupling fii. We are limited in 
our ability to do the former by our pulse timing resolu- 
tion, 5t, the uncertainty in pulse length Ti. 

Most qubits in the lattice will be far from the target 
qubit so that /(r, z) is small (Eq. ( 19 1), and thus will not 
experience any Stark shift due to the addressing beam. 
The |0) <-^ |2) and |1) ^> |3) transitions are detuned 
from these unshifted qubit transitions by Aac, while the 
|0) ^^ |1) transition is detuned by 2Aac. Since the prob- 
ability P of transition for any given atom is small, we 
can treat the transition amplitudes as independent, and 
calculate each transition probability using the Rabi for- 
mula, Eq. (20a), assuming that the coupling and pulse 

TTr 



time are chosen such that sin 



A2: 



= and 



"1 I " 2 

T[ ~ 7r/rii, where T[ ~ Ti/3 is the time required for one 
"leg" of the single qubit gate (see above and Figure |2| . 
Since we wish to minimize Aac for the purposes of reduc- 
ing other types of errors discussed below, we note that 
the smallest value of A for which the former condition is 
satisfied is A = v3^i- This results in the off-resonant 



transition probability estimates 
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With these estimates we can now ask, what is the corre- 
sponding EPG due to ofF-resonant transitions of all non- 
target atoms? If we can simultaneously address an entire 
row of riA = N^^^ atoms in a lattice of N atoms, the EPG 
is N'^/^Pmo- Unfortunately, this is challenging in even a 
modestly-sized lattice, as the intensity of the beam is 
inhomogeneous along the beam axis, on a length scale 
set by the Ray leigh length zr- By combining equations 



(191 and (20a), with A = VSf^i ( 



1- 



we esti- 



mate that the error for an atom at distance z from the 
beam waist along the beam axis would be approximately 
1 — Pmo{z) — ^Tirfr- For typical parameter values, this 
limits us to addressing just a few lattice sites with a sin- 
gle beam before the error becomes large. In fact for our 
example parameters, the Rayleigh length is zji — b.\ /im, 
and the error becomes of order unity at just one lattice 
site away from the beam waist. 

If we cannot simultaneously address entire rows with 
a single addressing beam, the effective EPG will scale as 

N ( I -^ j . Although such scaling of an error mech- 
anism would preclude scalable fault-tolerant quantum 
computation for an arbitrarily-large system, in practice it 
should not prove very restrictive for lattices with moder- 
ate numbers of qubits. For example, for microwave pulses 
of the appropriate frequency, 6t can be on the order of 
one cycle, or 10~^° s. For a single-qubit gate time of 
Ti = 10/J,s, this means a lattice of 10^ atoms could have 
an EPG of less than 10~^ due to ofF-resonant transitions. 
Scaling implications are discussed further in Section |III[ 



2. Addressing beam-induced heating 

The far ofF-resonant Gaussian addressing beam used 
to perform site-selective single qubit gates contributes 
harmonic trapping and anti-trapping potential terms for 
the |2) and |3) states, respectively, and also adds addi- 
tional state-dependent anharmonic terms to the poten- 
tial experienced by the atom. These anharmonic terms 
can generate entanglement between motional and inter- 
nal degrees of freedom, as well as "heating" the atom to 
higher motional states. Perturbation theory shows that 
the most significant undesirable effect is due to the differ- 
ence between the harmonic components of the trapping 
potentials experienced by atoms in the auxiliary states 
1 2) and |3), relative to those experienced in qubit states 
|0) and |1). 

The overlap between the vibrational ground state of 
the |0) level and the first even vibrationally-excited state 
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FIG. 3: A series of "snapshots" of the atomic spatial and 
internal state wavefunction for an atom in a ID trapping po- 
tential undergoing a microwave pulse-based single qubit n ro- 
tation gate, as simulated by QSIMS. The four relevant in- 
ternal atomic states are labeled |0) through |3), as in Fig. 
[2J and the curves represent the wavefunction for the atom's 
center-of-mass coordinate. Simulation parameters have been 
chosen to exaggerate the gate errors so that they are visible 
in this figure. In (a), we see the initial state of the atom is the 
internal state |0) {F — 3, nip ~ 0), with a Gaussian center- 
of-mass wavefunction. In (b) and (c), we see transitions to 
the auxiliary levels |2) and |3). In (d) and (e), vibrationally- 
excited states become visible, particularly for the component 
of the total wavefunction in the state |3). In (f), we see the 
final state with a significant portion of the wavefunction not 
in the desired |1) state, showing instead noticeable entangle- 
ment between spatial and internal degrees of freedom (e.g., 
the "bumpy" shape of the wavefunction component for the 
1 3) level). 
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where 



of the 1 2) level is approximately ^t 

LUr is the characteristic trapping frequency defined above 
and wq is the addressing beam waist. From this, we 
can calculate the probability of transition to the first 
even vibrationally-excited state, using the fact that uir ^ 
^T^i and Ti ~ Tr/fii: 
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Here we have simplified Eq. (21 1 by assuming that, since 
lOt is of the same order as fii, the sin (...) term will be 



of order unity. 

We have tested these perturbative predictions by sim- 
ulating the single-qubit tt gate with QSIMS as described 
above, with initial condition |0). The simulations show 
that the probability that the atom will not be in the 
motional ground state and desired qubit state after com- 
pletion of the IT gate is 1 x 10"^. This is consistent with 
what is expected from application of Eq. (21 1 with our 
simulation parameters. 

A small amount of heating does not in itself destroy 
the qubit state. However, as the vibrational energy of 
the atoms increases, the probability of other types of er- 
rors increases, and so the atoms will periodically need to 
be re-cooled to the ground state. Optical cooling can be 
done directly if the qubit state were first transferred to a 
different location, or could potentially be done with the 
qubit "in place" through a mechanism such as sympa- 
thetic cooling. Analysis of such re-cooling mechanisms is 
beyond the scope of this paper. 



3. Addressing beam-induced Raman Scattering 

We now determine the Raman scattering rate for 
the target atom during a single-qubit microwave gate 
(Rayleigh scattering at this wavelength is negligible) . Us- 
ing Eq. (pi, we can find the scattering rate in terms 
of E^. For e+ polarized lig ht at the magic wavelength 
Am, calculation using Eq. (A. 3 1 yields a polarizability 
of la^-U 



cording to: 



nm)| = 2.5 x lO^^C^ m^ J" - tor -""Us m 
the \F — 3, mp = 1) or \F = 4, nip = 1) state (note that 
the polarizabilities for these states have opposite signs 
and that the polarizabilities for the mp = states are 
essentially zero at Am)- Since the AC Stark shifts of the 
target atom \F = 3, m,p = 1) and \F = 4, mp = 1) lev- 

els are given by fiAac — — ^^(Am), we can express the 
scattering rate in terms of Aac as follows: 



T/h = 880 nm- 



coAq 



3.4 X 10"^ A„, 



nm) (22) 



We then obtain the corresponding Raman scattering er- 
ror per gate, Pmsi by multiplying the scattering rate T/h 
of Eq. ( 22 ) by the single-qubit gate time Ti . 
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There is also an effect due to the perturbation of 
the eigenstates of the |2) and |3) states. The matrix 
element between the unperturbed harmonic oscillator 
ground state and the perturbed first excited state is 

■5x 



^Sx = ^^^^\li^^^^- Assuming ^^^f^i < w^, the prob- 
ability of exciting to a higher motional state during the 
|0) ^^ |2) or |1) ^^ |3) transitions is then: 



raph 
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(24) 
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For typical parameter values, this second effect is of 
greater significance and we will neglect the former in com- 
parison with this. Note also that we have calculated the 
error only for one "leg" (i.e., transition) out of the three 
that compose the gate, and that the total error for the 
gate may be greater. 

We have simulated the one-qubit microwave pulse- 
based TT gate between |0) and |1) using QSIMS with an 
addressing beam position error S^ = 0.01 fim. We find 
that on completion of the gate, the probability that the 
atom will not be in the motional groun d st ate and desired 
qubit state 1 1) , is 2 x 10^^ . Using Eq. ( 24 ) with the same 



parameters as this simulation yields a value 1.3 x 10 ^. 
This is an estimate of the error in the |0) -^ |2) and 
|3) -^ |1) legs of the gate. For the |2) -^ |3) leg, we 
replace 6^ —^ 2Sx (because the perturbation due to the 
addressing beam has opposite sign for the |2) state versus 
the |3) state), to obtain an error of 7.8 x 10^^. Summing 
these three errors, we obtain an overall error estimate of 
7.8 X 10~^ for the complete gate, which is within a factor 
of 3 of the value obtained from the simulation. 



D. Two-qubit gates 



4- Addressing beam position error 

If the addressing beam is off-target by an amount S^ , 
Taylor expansion of (19) shows that the energy of the 



11) and 12) levels will be shifted by an amount 2hAac^- 
This decreases the |0) -^ |1) transition probability ac- 



In this section, we make a qualitative comparison of 
the two-qubit gate techniques most commonly mentioned 
in the literature, with an emphasis on analysis of their 
different implications for scalability. Most schemes either 
involve stationary qubits and long-range interactions, or 
movable qubits and short-range coUisional interactions 
between qubits. 
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1. Long-range interaction-based gates 

As with single qubit gates, our choice of a large lattice 
spacing allows for two-qubit operations to be performed 
in a site-specific manner. If the atoms are to remain sta- 
tionary, a long-range interaction is required to perform a 
two-qubit gate. Furthermore, the interaction must some- 
how be controllable. Dipole-dipole interactions between 
atoms excited into Rydberg states are a promising can- 
didate, and many variations on this idea have been pro- 
posed [lEIlI^IMlEIllMllMlCQJ. 

In one version of the Rydberg gate @], two nearby 
atoms are conditionally excited into Rydberg states via 
a coherent process [71 . If both atoms were initially in 
the |0) state, they are both excited into the Rydberg 
state where they interact via a dipole-dipole interaction, 
resulting in accumulation of a phase on the two-qubit 
state. They are then de-excited via another series of Ra- 
man pulses. Since the interactional phase accumulation 
occurs only in the case where both atoms are initially in 
the |0) state, this is effectively a CPHASE gate. Esti- 
mates of both the speed and the maximum possible fi- 
delity for Rydberg gates are reasonably promising |Hl51j. 
with some putting the error rates achievable on the order 
of 10^"^ to 10^^ j69j. However, because of the inherent 
long-range nature of the interactions in this gate, when 
running these gates between multiple pairs of qubits in 
parallel, it is essential to take into account the effects of 
cross-talk between different qubit pairs, i.e., the dipole- 
dipole interactions between qubits from different pairs. 
To make a rough estimate of the degree of paralleliza- 
tion possible in the presence of such cross-talk, consider 
a three-dimensional lattice of atoms. Suppose an exter- 
nal static electric field is applied to induce a "permanent" 
dipole moment in the Rydberg-state atoms [T: the level 
shift due to the resulting dipole-dipole interaction falls off 
as 1/R^. The fidelity error due to cross-talk will therefore 
scale roughly as l/R^, where R is the distance between 
different pairs of atoms that are simultaneously involved 
in Rydberg gate operations. We take the fidelity of a 
Rydberg gate performed between atoms one lattice spac- 
ing apart as unity, for reference. Then, taking a value 
of 10~^ for the fault tolerance threshold for gate errors 
(a value intermediate between different estimates of the 
threshold for the case of local gates [HI Hj), this impKes 
that the two-qubit gate cannot simultaneously be per- 
formed on multiple pairs of atoms within approximately 
ten lattice sites of each other. In a three-dimensional lat- 
tice, this geometric constraint limits us to simultaneously 
performing approximately one two-qubit gate per several 
hundred atoms. This in turn implies that the storage er- 
ror rate during the two-qubit gate time would have to 
be about two orders of magnitude lower than this fault 
tolerance threshold (i.e., ^ 10~* for the above example), 
to avoid accumulating additional idle time errors on the 
qubits not involved in the gates. 

It is possible to mitigate this cross-talk limitation by 
using an interaction with more limited range, such as 



the van der Waals interaction present between Rydberg 
atoms when there is no hybridizing static electric field 
[72] . In general, interaction strengths for the zero ex- 
ternal field case scale as 0{n* /R^) (as is typical for van 
der Waals-type interactions), although there are excep- 
tions [73] . This implies that cross talk errors would scale 
as 1/i?^^, which allows for roughly one simultaneous two- 
qubit gate per three dozen atoms. This imposes more 
modest constraints on storage error rates — the storage 
error rate per two qubit gate time would then only need 
to be about one order of magnitude below the threshold 
value of 10^6. 



2. Collisional gates 

Another method of avoiding the limitations of long- 
range interactions is to bring pairs of atoms close together 
and use short-range interactions or collisions to imple- 
ment two-qubit gates. A variety of such gates have been 
proposed and analyzed [51|71|751|7Sl[771[7Hl[7i[HnilHI]- 

We note that, in the context of the system discussed in 
this paper, it would first be necessary to transfer atoms 
from the qubit states |0) = |F = 3, mj;- = 0) , |1) = 
\F = A,mF=Q) to states that experience a state- 
selective trapping potential (i.e., states with non-zero 
mp), so as to allow pairs of atoms to be translated 
towards each other and brought together. For exam- 
ple, for lattice light with A = 800 nm and a par- 
ticular circular polarization, the states \F — i^mp — 1) 
and \F = 3, nip — —1) have polarizabilities that differ by 
more than 8 %. Specific atoms can be transferred into 
these states by microwave or Raman transitions, analo- 
gous to the single qubit gates discussed earlier. Alter- 
natively, an entire plane of atoms can be transferred si- 
multaneously using the microwave pulse method if the 
addressing beam is replaced by an inhomogeneous mag- 
netic field. Once the atoms are in the appropriate states, 
the atoms can be selectively moved by changing the rela- 
tive polarization of one of the lattice beams, as described 
by Vala et al [S2] and Weiss et aZ |7] , allowing the atoms 
to be brought together to perform a gate |5]- 

Due to the necessity of physically moving atoms 
around the lattice, collisional gates are likely to be much 
slower than long-range interaction-based gates. Even in 
the case of a "fast approach" , where the translation of 
the atoms is not adiabatic, gate times are still one to two 
orders of magnitude slower than the characteristic trap 
period of the lattice site's potential well iTZ]. Estimates 
of the maximum fidelity possible with collisional gates 
also tend to be lower than maximum fidelity estimates 
for Rydberg gates [7i] . 

Despite these drawbacks, collisional gates appear more 
suited to large-scale quantum computation than long- 
range interaction-based gates. Collisional gates can eas- 
ily be performed on a massively parallel scale, particu- 
larly when used in the context of cluster-state quantum 
computing [H5J [55] . In cluster-state (also known as "one- 
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way") quantum computing, two qubit gates are used only 
in the initial preparation of a large entangled "cluster" 
state. The computation itself is then effected via single- 
qubit measurements in a variety of bases, or equivalently, 
single qubit rotations followed by measurement in a par- 
ticular basis. One recent scheme for cluster-state quan- 
tum computing offers fault tolerance thresholds as high 
as lO^'^ for local depolarizing errors and lO^'^ if there are 
also errors in preparation, gates, storage, and measure- 
ment [H]. One can easily imagine building such a cluster 
state with atoms in a three-dimensional lattice by per- 
forming coUisional gates in parallel on entire planes of 
atoms at a time. 



III. ANALYSIS 
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FIG. 4: (Color online) A 3D plot of gate error (e) versus 
lattice spacing (a) and beam waist (wo) for the Raman single 
qubit gate, with lattice depth Ul = 500 fiK and detuning 
lAil =27r-5 THz. 



Having quantified the various sources of error in the 
previous section (summarized in Table m , we now make 
use of these results to analyze the extent to which single- 
qubit operations can be performed in parallel in an op- 
tical lattice. We shall first determine an estimate of the 
lower bound on the error rate due to single qubit gate 
errors for the case of a blue-detuned lattice with Raman 
transition-based single qubit gates. We then discuss the 
implications of these gate errors and lattice based storage 
errors, together with lattice addressability issues for the 
scalability of quantum computation in three- and two- 
dimensional lattices. Following this analysis of Raman 
gates, we then discuss how the corresponding arguments 
apply in the case of the microwave single qubit gate. 



A. Raman-based single qubit gate errors 

We sum the EPG for each relevant error mecha- 
nism to obtain a total error rate for the gate. Using 
(l7| and the expression for the intensity of a Gaussian 
beam of power Pr, we see that we can replace \^r\ by 
(8.3-1012 W^i m^ s-2)^;^|^L. We also require that the 

Raman beam laser power Pn not exceed some realistic 
value Pmax- We then numerically minimize the resulting 
expression for the total Raman gate error over the param- 
eters a. Wo, Ai, Ul, and flu (see Raman gate terms in 
Table [11 Carrying out this procedure, we find that even 
if we choose an exceedingly high value of Raman laser 
power, Pmax = 10 W and allow the other parameters to 
attain unrealistic values, we do not obtain a minimum 
EPG below 10"^ for the combined total Raman gate er- 
ror mechanisms listed in Table [ij In fact , this minimum 
value is already unrealistically low since at these laser 
powers ionization of the Cs atoms would also be expected 
to come into play. 

If we now use a more realistic value of Raman laser 
power of Pmax = 10 mW, and impose realistic constraints 
on the other parameters, corresponding to a < 10 /im, 
Ul < 500 ^K, I All < 27r • 5 THz, we obtain a minimum 



EPG of 1 X 10~^. This minimum occurs when the afore- 
mentioned parameters reach their constraint values and 
wq = 5.0 /im. This EPG is roughly an order of mag- 
nitude better than can be achieved with our reference 
parameters of a = 5 um, Ul = 200 /iK from Section |T] 
A plot of the EPG as a function of a and wq is shown in 
Fig. |4J In this Raman-based single qubit gate scenario, 
polarization effects are the dominant source of gate er- 
rors, and the minimum achievable EPG is consequently 
most sensitive to changes in the lattice spacing, a. A 
doubling of Cmax to 20 /xm (and appropriate adjustment 
of Wq) results in an EPG that is approximately threefold 
smaller. It is also worth noting that although the EPG 
is sensitive to the laser power at the unrealistically high 
values discussed above, it is not sensitive to Pmax in this 
regime. 

This analysis suggests that one consider whether it is 
possible to further increase the lattice spacing a beyond 
the maximum value 10 fim specified above. While in 
principle, this would improve (reduce) the EPG, we note 
that the laser power at a particular frequency required 
to generate an optical lattice of given depth and a given 
number of atoms scales with the square of a. Thus scaling 
to larger lattice sizes will necessarily entail significant in- 
creases in laser power that may not be realistic. We note 
that while larger lattice spacings are beneficial for the 
single qubit Raman gate, they will result in slower gate 
times for two-qubit gates, as a result of the weaker in- 
teraction strengths for long-range interaction-based gates 
such as the Rydberg gate and of longer travel distances 
in the case of coUisional gates. 



B. Lattice size and scaling with Raman- based 
single qubit gates 

We now examine the implications of this analysis of 
the total Raman-based single qubit gate EPG for lattice 
size and scaling. 
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Source 



Section 



EPG 



Raman scattering (blue-detuned lattice) 


II A 1 




Raman scattering (red-detuned lattice) 


II A 1 


ri;^i^f/4-.,i?§)^|^ 


Neighbor atom errors (Raman gate): Pm 


II Bl 




Spontaneous emission (Raman gate): P^a 


II B 2 


7r 


2|Ai|t 


AC Stark Shifts (Raman gate): P^a 


II B 3 


4 h'-a^ /'e-lQlO s-l\^ 
371^ mwf.UL \ ^1 ) 


Atomic motion-reduced pulse area (Raman gate): Prpa 


II B 4 


1 r2^2 

12 mw^Ui^ 
r, no / 27r \ h rrl/2 


Detuning Doppler shift (Raman gate): Prm 


II B 4 


^■■^^\x^n^) m3/2„t^L 


Polarization effects (Raman gate): P^p 


II B 5 


nxla 1 


T^^O ^^■2mUL 


Laser intensity noise (Raman gate): Pri 


II B 6 


C^ fi 1 n^'^ m^ "^^"1 1 


^ • ' |Ai|tu^ 



TABLE I: Error Per Gate (EPG) due to various sources. It is assumed that atoms are in the vibrational ground state (i.e., 
rix = fly = G). The first two sources listed produce storage errors, while the rest cause gate errors. 



Sectii 



EPG 



Numerical EPG 



Raman scattering (bluc-dctuncd lattice) 


II A 1 


Ti 


N 
^A 


.cso Ur,(u:^,E^) ,( 


-i) 


2.4 X 10"'' 


^^L V =- k( 


"i)| 


Off-resonant transitions (microwave gate): Pino 


Pmh 


II CI 






1 ^'^l.rn.<> 




4.3 X 10"** 


Addressing beam-induced heating (microwave gate) 


lie 2 


1 X lO"** (*) 


64,4 Tful^fj 

3.4 X 10-" AacTi 

VTsaAL.52„5„3/2 
,3 2„5/2 4 


Raman scattering (microwave gate): P^s 


lie 3 


5.2 X 10"^ 


Addressing beam position (microwave gate): P^p/j, 


lie 4 


2 X 10"^ (*) 



TABLE II: Error Per Gate (EPG) due to various sources. In calculating a numerical value of the EPG, we used the parameters 
A^ — 10^, Ha = WO, Sti = 10^^" s,5x = 0.01 fira,LJL ~ 27rc/800 nm, with the other parameters as defined as in Sec. II C 
(a = 5 /im, [/l — 200 /iK, Ti — 76 ^s, Aac ~ 0.2 MHz,uio = 1.2 fim) Where indicated above by (*), we use numerical values 
obtained from simulations rather than from the analytic approximations listed. We treat heating "events" as errors, although 
it is worth noting that if a re-cooling mechanism is implemented, the effect of such errors can be reduced. 



3D lattices 



With three sets of 10 W laser beams, we can produce 
a 3D lattice of ~ 100 x 100 x 100 sites, a lattice spacing 
of 10 ^m, and a lattice depth of ~ 500 fj.K, by tuning the 
laser to ~ 851.7 nm (which is very close to the D2 tran- 
sition at 852.1 nm). With this small a detuning, there 
is a substantial and undesirable ~ 3% difference in trap 
depth for atoms in the |0) state versus the |1) state. This 
mismatch could cause entanglement of motional and in- 
ternal qubit degrees of freedom and should be avoided. 
Rayleigh scattering would also cause rapid heating of the 
trapped atoms, requiring frequent re-cooling operations. 
There is thus little if any room to decrease detuning fur- 
ther. Consequently, larger lattices with the same trap 
depth that might allow a lowering of the EPG according 
to our analysis above would be possible only by making 
substantial increases in laser power. 

What are the implications of these results for the com- 
puting power of a 3D optical lattice-based scheme? Ob- 
viously, given some limit on the available lattice laser 
power, the lattice size and thus the number of physical 



qubits will be limited, with the number of qubits scaling 

3 /2 

as ~ Pmax. Another constraint is the degree of paralleliz- 
ability of performing gates in an optical lattice. As dis- 
cussed earlier in Section [ll] fault-tolerance threshold the- 
orems assume maximal parallelizability |13j . which im- 
plies that all or nearly all qubits need be addressable 
simultaneously (ha — N, where ha is the number of ad- 
dressable qubits and N is the total number of physical 
qubits). In the case of 3D addressable optical lattices, 
one might imagine using micro-lens arrays |84| to allow 
simultaneous focusing of many addressing beams. It may 
initially appear possible to address ~ 7V2/3 g^^Qj-^g jn a 3D 
lattice with this approach, using arrays of ~ N^^"^ x A^^/^ 
lenses placed adjacent to the lattice. However, these 
lenses must be able to focus on sites deep inside the lat- 
tice. Since the linear size of the lattice and thus the 
distance to such sites scales as N^/^ , the diameter of the 
lenses must also scale linearly with N^/^. This implies 
that the number of addressable atoms nA in a 3D lattice 
is more or less constant, and does not scale with N . We 
note that 2D optical lattices (see below) do not neces- 
sarily have the same limitation, since the distance from 
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3D lattice and not a 2D architecture. The use of a 3D 
lattice implies a somewhat better threshold due to the 
reduced average distance between qubits in 3D. Overall, 
this analysis suggests that Raman gates may indeed meet 
the fault tolerance threshold within the assumptions of 
our analysis (i.e., realistic laser power), if parallelizability 
limitations can be overcome. 



FIG. 5; (Color online) A method for addressing some or all 
of the atoms in a single 2D plane of a 3D optical lattice. The 
circles on the right represent atoms in an optical lattice (only 
a cross-sectional plane of atoms is shown). The stars on the 
left represent "point sources" of light, which could consist of 
lenses on the end of individual optical fibers, with each fiber 
coupled to a light source and independently controlled. The 
lens in the center focuses the light from the point sources 
onto corresponding atoms in a plane (i.e., the vertical plane 
perpendicular to the page surface). By adjusting either the 
location of the lens or the array of point sources, different 
planes of atoms could be targeted. 



the micro-lens array to the focal plane would be constant 
rather than scaling with the linear lattice dimension. 

Another approach for 3D lattices would allow address- 
ing of up to N^/^ sites in parallel, using a single large lens 
and a 2D array of light sources (see Fig. p|. It is in princi- 
ple possible to address up to all ~ N atoms using a setup 
similar to Fig. [Slbut with the "point sources" replaced by 
a spatial light modulator. However, the resolution (total 
number of pixels) of the spatial light modulator would 
have to scale with N, posing considerable challenges for 
implementing such an approach. 

Given these constraints, how large a computation 
could be performed with a 3D lattice using current tech- 
nology? To address this question, let us consider the 
aforementioned 100 x 100 x 100-site lattice (generated 
with three 10 W beams). In this scenario, the per-atom 
Raman scattering rate for lattice light is approximately 
0.4 s~^. If we use 1 W Raman beams, we can perform 
single qubit gates in approximately half a nanosecond on 
up to 10^ qubits in parallel. This means that the single 
qubit lifetime due to Raman scattering of lattice light 
is approximately 10^ time-steps, and that 10^ gates per 
qubit could be performed in that time. 

We see that in such a scenario, the storage error rate 
due to Raman scattering rate of lattice light are almost 
two orders of magnitude smaller than the gate error rate, 
which we calculated above to be EPG = 1 x 10~^. It 
is useful to compare these numbers with the most opti- 
mistic rigorous threshold result for local gates in a 2D 
architecture, which is currently 1.85 x 10~^, with storage 
errors an order of magnitude smaller |12j . This compar- 
ison is complicated by the fact that we are considering a 



2. 2D lattices 



2D lattices are subject to the same basic size constraint 
in terms of laser power. This reduces the number of 
qubits accessible, relative to a 3D lattice formed with 
equivalent power. Thus, a 2D lattice with ^ 10® phys- 
ical qubits can be formed in a 1000 x 1000-site lattice 
with spacing a ~ 10 /im and depth Ul ^ 500 /iK. It is 
straightforward to show from the analysis above that this 
does however now require kilowatts of power, even with 
the small lattice light detuning given earlier. A 100 x 100 
lattice with the same depth and spacing could be created 
with more realistic ~ 10 W laser beams, and this would 
then accommodate only ~ 10* physical qubits, in con- 
trast to the 10® that were possible with the 3D lattice 
above. 

2D lattices are more attractive from the perspective of 
simultaneous addressability. A micro-lens array or the 
method described in Fig. |5] can in principle be used for 
addressing lattice sites in a 2D lattice in such a way as 
to allow the number of simultaneously-addressable sites 
UA to scale linearly with the total number of of lattice 
sites N . This is possible for a 2D lattice because such 
methods allow for parallel addressing of at most all the 
atoms in one plane, and a 2D lattice of course consists 
only of a single plane of atoms. It was recently shown 
that multiple standing-wave fields can also be used to 
address a periodic subset of the atoms |33] . 

Gates performed in a 2D lattice will also not be subject 
to all the same error mechanisms that we would find in 
a 3D lattice. In particular, in a 2D lattice it may be 
possible to use a first-order Doppler-free configuration 
for the Raman beams, which would substantially reduce 
the detuning Doppler shift error. Nevertheless, for the 
Raman-based single qubit gates, polarization effects are 
the dominant error mechanism, and these remain more 
or less unchanged by the switch to a 2D lattice. Thus 
we still obtain an EPG value of 1 x 10~^ using the same 
parameters given above for the 3D lattice analysis. As 
noted there, if no restrictions on laser power existed, we 
could increase the lattice spacing a and the beam waist 
Wo to reduce the EPG. In practice, laser power is likely 
to be the limiting factor on scaling of 2D lattices of the 
type considered here. 
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C. Lattice scaling with microwave-based single 
qubit gates 

Our scaling analysis above has focused on the Raman 
single qubit gate. The microwave gate is subject to the 
same fundamental constraint on scaling in a 3D lattice, 
since it also requires a focused beam to address indi- 
vidual lattice sites. Microwave-based gates are typically 
substantially slower than Raman gates, taking on the or- 
der of tens of microseconds. This implies that lattice 
light scattering rates would need to be lower for a fault 
tolerant quantum computer employing microwave-based 
gates, which could be achieved with larger lattice light 
detuning and increased laser power. On the other hand, 
lattice spacings need not be quite as large for microwave- 
based gates, reducing the laser power required for a lat- 
tice with a given number of sites. As the former effect is 
larger than the latter, we expect systems of a given num- 
ber of qubits employing microwave-based single qubit 
gates to require more power than systems of a similar 
number of qubits using Raman single qubit gates. 

For example, using the parameter values given in Table 
|TT]for a microwave single qubit gate-based scheme, three 
75 W lasers would be required to create a 100 x 100 x 100 
lattice of depth Ul = 200 /iK with spacing a = 5 /xm and 
a lattice laser wavelength of Xl — 800 nm; despite the 
large detuning, atoms in such a lattice would still suffer 
from a large effective EPG of 2.4 x 10^^ from Raman 
scattering of lattice light, due to the slow gate times and 
limited parallelizability of microwave gates. Note that 
this error rate for scattering of lattice light is substan- 
tially larger than the total error rate of ~ 7 x 10~^ for 
effects intrinsic to the microwave gate itself. By com- 
parison, a Raman single qubit gate-based setup with the 
same number of lattice sites could make due with laser 
powers of 10 W per laser and still have an effective EPG 
for lattice light-induced Raman scattering that is five or- 
ders of magnitude smaller than that for the microwave 
gate scenario (see Sec. Ill B 1 for details) . 



While Raman single qubit gates are likely to be prefer- 
able for scaling to very large system sizes, microwave- 
based gates may nevertheless prove easier to implement, 
due to the simpler optical requirements. In particular, 
beam alignment is simplified in the case of microwave 
gates, and since they require only a single laser, they 
will be less sensitive to alignment errors than the Raman 
single qubit gates. 



IV. CONCLUSION 

We have investigated the effects of single qubit gate 
errors and lattice storage errors for neutral atom qubits 
trapped in addressable optical lattices and analyzed the 
role of these errors and of other technical factors such as 
finite laser power in limiting the size of a quantum com- 
putation. We find that under realistic current limitations 
on laser power, computation with lattices containing up 



to 10^ qubits in three dimensions and up to 10^ qubits 
in two dimensions may be achievable. Considerations 
of parallelizability and optical access are seen to impose 
additional limits on the scalability of such quantum com- 
putation with individually addressed gates. These con- 
straints are more severe for a 3D lattice than for a 2D 
lattice, where it is comparatively easier to develop tech- 
nologies for addressing all or nearly all atoms in parallel 
(i.e., 71 A ^ N). 

Our quantitative scaling analysis has not included the 
effects of two qubit gate errors and measurement errors. 
For the single qubit gates, we compared the accuracy of 
implementation for two candidate gates: the stimulated 
Raman two-photon gate, and the microwave gate with 
AC Stark Shift addressing beam. 

We find that Raman-based single qubit gates can be 
implemented for ^'^^Cs in times on the order of a nanosec- 
ond with an error rate of ~ 1 x 10~^ and that microwave- 
based single qubit gates can be implemented in times on 
the order of 100 /iS with an error rate of ^ 7x10^^^. (Note 
that neither error rate includes scattering of lattice light, 
which is a more serious limitation for the microwave gates 
due to their slower gate times.) The microwave gates 
are simpler to implement (e.g., require less laser align- 
ment), but arc more severely limited by constraints on 
site-specific parallelization due to the global nature of 
the microwave pulse. Consequently, microwave gates ap- 
pear to be a viable intermediate option for testing single 
qubit gates and realizing small scale quantum algorithms, 
although Raman-based gates may be preferable for even- 
tual operation at the full capacity of addressable optical 
lattices. The gate error for Raman-based single qubit 
operations is very close to the current best estimate for 
the fault tolerant threshold for computation using local 
gates in a 2D architecture. The latter would be expected 
to be somewhat higher in a 3D architecture, suggesting 
that further innovations in error correction protocols may 
make reasonably large computations, i.e., with up to 10^ 
qubits, achievable with 3D addressable optical lattices. 

The focus of this work was a realistic analysis of tech- 
nological limitations on scalability of fault tolerant quan- 
tum computation with neutral atoms in addressable opti- 
cal lattices. Our analysis identified laser power and par- 
allel addressability as primary factors that will eventually 
restrict the number of qubits for this implementation of 
quantum computation. We hope that this detailed anal- 
ysis will stimulate similar investigations of physical and 
technological limits to other proposed implementations of 
fault tolerant quantum computation, each of which has 
very different limiting physical features whose effect on 
scalability needs to be examined in detail. 

It may also be possible to overcome the limits identified 
here for scalable quantum computation in 3D lattices, 
by making use of methods other than single-site focused 
beams in order to address large numbers of atoms in par- 
allel. Thus, entire planes of atoms can be addressed using 
a magnetic field gradient [5] . Alternatively, periodically- 
spaced arrays of atoms can be addressed with a sec- 
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ond optical lattice of larger lattice constant (a "super- 
lattice" ) [ISl ES] • Both of these techniques offer high de- 
grees of parallelization at the cost of flexibility in choos- 
ing which atoms are simultaneously addressed. Research 
on quantum cellular automata i86j and global control 
schemes [87] has shown that such a constraint on simul- 
taneous addressing does not necessarily preclude efficient 
quantum computation, but further work is needed to un- 
derstand the effect on error correction protocols and to 
determine fault tolerance thresholds for such schemes. 

Finally, we note that despite the limitations estab- 
lished here for performing large scale quantum com- 
putations, optical lattice-based schemes may also be 
very useful for simulating other quantum systems [5H] . 
Quantum simulations have different requirements than 
fault-tolerant quantum computation, generally requiring 
less stringent accuracy of quantum operations. Optical 
lattice-based quantum computers may also be useful for 
small-scale quantum information processing, such as a 
quantum repeater in a quantum key distribution net- 
work |89j . The long coherence time of optically-trapped 
neutral atoms is particularly valuable for this application. 
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magic wavelength, Am- We find it valuable to review 
proper methods for these calculations. 

For a quantum treatment of scattering, the relevant 
starting point is the Kramers-Heisenberg formula (A.ll, 



which can be derived using second-order time-dependent 
perturbation theory: 



dcTab = 
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where a is the fine structure constant, the indices a, b, 
and i denote initial, final, and intermediate states, x^a 
is the dipole matrix element and w^ — ^i — ^a the fre- 
quency for the a — > i transition, lo and uj' are the fre- 
quencies of the incoming and scattered photons, and e 
and e' are the polarization vectors of the incoming and 
scattered photons. In the most general case, the sum over 
discrete intermediate states i should be extended to in- 
clude an integration over continuum intermediate states, 
but, for our purposes we can safely neglect continuum 
states and other far-detuned intermediate states. Inte- 
grating over dU, and multiplying by 2 to account for the 
two possible independent photon polarizations, we ob- 
tain the cross-section for a state a to scatter into a state 



h (A.2): 
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Raman scattering cross section is given by 
CTraman.a = S&^a '^ab' whilc thc Raylcigh Scattering cross 

section is CTrayloigh.a = <7aa- 

Dynamic polarizability for an atom in state a (e.g., for 
an AC Stark Shift) is given by (A.3), which can be de- 



rived using second-order time-independent perturbation 
theory: 



APPENDIX: SCATTERING CROSS SECTIONS 

AND DYNAMIC POLARIZABILITIES IN 

GROUND-STATE ALKALI ATOMS 

In analyzing optical lattice-based quantum computing, 
one frequently has to calculate optical-frequency scatter- 
ing cross sections or polarizabilities for an alkali atom in 
a particular hyperfine sub- level of the ground state (e.g., 
in ^•^^Cs, 6^Si/2, F — 'd, mp — 1). Such calculations 
can be performed using well-known formulae; unfortu- 
nately, due to thc tedious and somewhat subtle nature 
of these calculations, various approximations are some- 
times employed as shortcuts. Such approximations are 
valid far from resonance, but result in significant errors 
when close enough to resonance to resolve the hyperfine 
structure, as is the case for an addressing beam at the 
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For both (A.2 1 and (A.3 1, the matrix elements we need 
to calculate are between a hyperfine sub- level \F mp) of 
the Si/2 ground state, and a hyperfine sub-level \F' m'p) 
of either the P1/2 or P3/2 (-Di-line or Z?2-line) excited 
states. Such matrix elements can be calculated using 
thc Wigner-Eckart Theorem and measured values for the 
appropriate reduced matrix elements |60l I90|. The ex- 
pression is given in the following equation (A.4|: 



(A.4) 



{F mp\^ ■ e*^\F' m'p) 

= (J||x||J')(-l)^+^+™^\/(2F-hl)(2F' + l)(2J-hl) 
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where eo = z and e±i = =F (x ± iy) /V^ (note that e±i 
corresponds to cr^-polarized Ught, and also that e±i = 
—e*i). Care must be taken with polarization vectors, as 



in general, x^a • e ^ Xai ■ £■ Thus, a particular intermediate 
state i may contribute to the resonant term but not the 



anti-resonant term of (A. 2 1 or (A.3), or vice versa 
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